Nonlinear c-axis transport in Bi 2 Sr 2 CaCu208+i from two-barrier tunneling 
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Motivated by the peculiar features observed through intrinsic tunneling spectroscopy of 
Bi2Sr2CaCu208+i mesas in the normal state, we have extended the normal state two-barrier model 
for the c-axis transport [M. Giura et al., Phys. Rev. B 68, 134505 (2003)] to the analysis of dl/dV 
curves. We have found that the purely normal-state model reproduces all the following experimental 
features: (a) the parabolic ^-dependence of dl/dV in the high-T region (above the conventional 
pseudogap temperature), (b) the emergence and the nearly voltage-independent position of the 
"humps" from this parabolic behavior lowering the temperature, and (c) the crossing of the abso- 
lute di /dV curves at a characteristic voltage V x . Our findings indicate that conventional tunneling 
can be at the origin of most of the uncommon features of the c axis transport in Bi2Sr2CaCu20s+«. 
We have compared our calculations to experimental data taken in severely underdoped and slightly 
underdoped Bi2Sr2CaCu208+5 small mesas. We have found good agreement between the data and 
the calculations, without any shift of the calculated dl/dV on the vertical scale. In particular, in 
the normal state (above T*) simple tunneling reproduces the experimental di /dV quantitatively. 
Below T* quantitative discrepancies are limited to a simple rescaling of the voltage in the theoretical 
curves by a factor ~2. The need for such modifications remains an open question, that might be 
connected to a change of the charge of a fraction of the carriers across the pseudogap opening. 



PACS numbers: 74.25.Dw,74.25.Fy,72.80.-r,74.72.Hs 

I. INTRODUCTION 

Tunneling is an inherent phenomenon in highly 
anisotropic cuprate superconductors— Due to their 
layered structure of conducting sheets intercalated by 
insulating blocks, superconductivity in these compounds 
is often described in terms of thin superconducting layers 
coupled via Josephson tunneling. Similarly, tunneling 
between the layers is expected to affect the transverse 
(c-axis) conductivity in the normal state, in competition 
with thermally activated interlayer hopping. If the 
interlayer barriers are not dramatically high, the latter 
process is a natural consequence at sufficiently high 
temperatures. 

In this sense, measurements of the c-axis conductivity 
add to the various other probes of the transverse 
quasiparticle transport: intrinsic tunneling spectroscopy, 
scanning tunneling spectroscopy, angle-resolved photo 
emission spectroscopy (ARPES). In particular, when an 
intrinsic tunnel junction exists, p c is a good measure of 
the density of states (DOS) around the Fermi level (low 
energy, in the spectroscopic language). 
Experimental investigations of the transverse properties, 
such as the differential conductivity dl/dV and ARPES 
(which probes the c-axis properties at selected momen- 
tum), reveal the existence of different energy scales and 
their evolution with doping and temperature. The most 
evident finding is the appearance of sharp features at the 
superconducting transition ("peaks") at typical energies 



As. The temperature evolution of Ag which, starting 
from zero, increases upon decreasing T from T c , ensures 
the identification with the superconducting gap. 
Much less obvious is the appearance of weaker features 
("humps") at larger energies A M , which is often identi- 
fied with the pseudogap. In this case the temperature 
dependence is much less pronounced. Sometimes a 
constant A pg (T) is observed decreasing T, while only 
the number of carriers is seen to decrease. The interpre- 
tation of the connection between Ag(T) and A pg (T) is 
still a very debated subject^ and a crucial issue for the 
understanding of high-temperature superconductivity. 
Measurements of the differential conductivity dl/dV 
gave a substantial contribution in bringing the mentioned 
features into a clearer light, at least on experimental 
grounds. The superconducting gap has been identified 
with the voltage V p where the peaks in dl/dV appear. 3,4 
A very weak temperature dependence of the bias voltage 
where the humps appeared, Vh, has been consistently 
reported in several papers4^^ In addition, it was 
shown that Vh was insensitive to a magnetic field, 
while V p was reduced to zero by large magnetic fields, 8 
suggesting a different nature of the phenomena involved. 
For what concerns the description of the data, the usual 
procedure has been to include some specific shape of 
the DOS in the model to fit the shape of the dl/dV 
curves. In particular, below T c it is now customary to 
use some BCS-like expression of the superconducting 
gap to obtain estimates of the maximum gap and of 
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the scattering rate from fitting of the dl/dV curves^ 
This approach has been extended to the fit of the 
humps. Taking two different gaps the dl/dV at 4 K 
were reproduced. 14 A continuation of the fitting of the 
dl/dV above T c with a nearly T— independent gap has 
been also reported.— The T— independence of Vh (Ref0) 
and, consequently, of the pseudogap if one identifies 
the experimental feature with a gap in the DOS, would 
imply that A pg is an energy scale for a crossover instead 
that for a true transition. 

Proceeding toward the normal state, we note that 
another interesting feature appears. The dl/dV curves 
at high T (and sufficiently large doping) present a very 
anomalous downward parabolic behavior fi&L which has 
been treated up to now as an unspecified "background" . 
In this paper we address the issue of the normal state of 
layered superconductors, focussing on P^S^CaC^Os+s 
(BSCCO), where most of the experiments have been per- 
formed. We build up a model for the c-axis dl/dV based 
on the presence of two different energy barriers, where 
each one gives rise to tunnel and thermally activated 
currents i 16 ! 17 ' 18 We show that this simple model, based 
solely on the existence of two energy barriers and with- 
out any particular features in the DOS, can reproduce 
many aspects of the experimental dl/dV. In particular, 
the parabolic behavior is reproduced quantitatively. The 
presence (in overdoped BSCCO) or absence (in under- 
doped BSCCO) of such high-T parabolic dl/dV finds an 
explanation in the simulations with the different balance 
between tunneling and thermal activation, depending 
on the height of the energy barriers. Moreover, the 
balance between thermal and tunnel currents across the 
barriers gives rise to humps in the dl/dV, suggesting 
that the humps arise from purely normal-state tunneling 
phenomena. The nearly constant position of the humps 
as a function of the temperature in experimental data, 
as well as the crossing of the dl/dV at a characteristic 
voltage in underdoped samples, are reproduced by the 
simulations. Finally, a simple rescaling of the voltage 
scale by a factor ^2 is sufficient to bring the simulations 
in accurate quantitative agreement with experimental 
data taken in BSCCO mesas. 

The paper is organized as follows. In Section [II] we recall 
the two-barrier model for the c-axis conduction, fo- 
cussing on the extension of the model to the description 
of nonlinear dl/dV vs. V curves. We find that several 
experimental features are recovered by the simple model. 
In Section IIHI we discuss data taken in BSCCO mesas 
in terms of the two-barrier model. A short summary is 
presented in Section HVl 



II. dl/dV IN THE TWO-BARRIER MODEL 

It is generally believed that the c-axis conduction 
in BSCCO is determined by the layered structure. In 
particular, tunneling certainly plays a major role. To ad- 



dress this issue we have introduced and develope d 16 ' 1 ^ 18 
a phcnomcnological model for the out-of-plane resistivity 
p c in layered superconductors. In particular, we have 
proposed that the two-layer structure of BSCCO could 
give rise to two barriers for the c axis conductivity. The 
overall conductivity would then be the composition of 
thermal and tunnel transport across each of the barriers. 
A sketch of the barriers considered is presented in Figure 
[TJ Clearly, this model directly applies to the true normal 
state, i.e. the state where no other energy scales appear 
(due to, e.g., the opening of the pseudogap or of the 
superconducting gap 1 ^). More specifically, no particular 
ad hoc shape of the density of states is considered. 
Regarding the scattering process involved, we have 
considered that tunneling across the barriers could be 
reduced due to incoherent in-plane scattering, along 
the lines of Ref.20. This process reduces the average 
interlayer tunneling hopping rate, thus leading to an 
increased charge localization in the (a, b) planes. In 
particular, in-plane inelastic phononic scattering given 
by the scattering time r reduces the tunneling matrix 
element t c to (2r/h)t 2 <C t c . In the linear response 
(vanishing voltage), the resistivity in the two-barrier 
model is then given by the simple expression ; 16 i 18 
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where the first and second term in square brackets repre- 
sent tunnel and thermal activation over the ith barrier, 
respectively. Thus, the model represents a series of two 
elements, each of them is a parallel of two channels. Scat- 
tering processes enter through r _1 cx ap a b.n (in- plane 
scattering) and b (out-of-plane scattering), where p a b.n 
is the (a, b) plane resistivity in the normal state, d\ = 3 
A and g?2 = 12 A are the spacings between the Cu02 lay- 
ers, and the height Aj of the i-th barrier enters through 
+2 1 - 
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For small applied voltages one can take a rectangular 
shape for the energy barriers without loss of generality, 
as depicted in Figure [1] upper panel. With this simple 
model we were able to reproduce the doping dependence 
of p c {T) in BSCCO quantitatively, and in particular the 
evolution from a nearly linear p c {T) in overdoped samples 
to the "semiconducting" behavior in underdoped sam- 
ples. In particular, the model was able to reproduce the 
local minimum of p c (T) which appears at optimal dop- 
ing and slight underdoping,- 1 — The details, including the 
crossover from thermal to tunnel transport by varying 
the temperature and the height of the barrier, have been 
worked out and discussed previousl y 16 ' 18 In the follow- 
ing we focus on the extension of this model to nonlinear 
transport. 

The problem of nonlinear tunneling is a long standing 
issue in solid state physicS ) 22 i 23 ' 24 i 25 ' 26 i 27 and it is still 
far from being fully solved^ Due to the extreme com- 
plexity of the problem, we use an approach as simple as 
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possible to incorporate the effect of finite voltages into 
our model. In particular, we have to incorporate (a) 
the deformation of the energy barriers due to increas- 
ing voltage, (b) the existence of two barriers "in series" , 
which implies noncqual, V^-dependent voltage drops on 
each barrier, leaving unaffected (c) the role of the in- 
plane scattering, which is a physical foundation of the 
linear model employed up to now. 

In a nonlinear problem, it is expected that the ap- 
plied voltage affects the shape and height of the energy 
barrier^ This effect cannot be taken into account by 
rectangular barriers. Among the many possible choices, 
we have found that parabolic shapes for the energy bar- 
riers led to an analytical, albeit complicated, expression 
for the tunneling current. We checked that the results 
did not change appreciably with other barrier shapes, 
such as trapezoidal or quartic or higher order polynomi- 
als. While the simulations were not very sensitive to the 
detailed shape of the barriers, two requirements emerged 
from extensive simulations: first, the barriers needed to 
vary their height with the applied electric field, that is 
rectangular barriers were found to be unsuitable for an 
even qualitative reproduction of the experiments. Sec- 
ond, we found that the existence of two different bar- 
riers ensured a nearly quantitative reproduction of the 
experimental dl/dV. Thus, we chose the energy profile 
sketched in Figure [1] (middle panel) to represent the two- 
barrier nonlinear model, with: Ui{x) = U C i — ^(x — ^j) 2 in 
the interval < x < di, and zero otherwise. With respect 
to the previously defined Aj, one has J7 C> , = A, + Ep and 
a = 8U C i/df. It is important to stress that the choice of 
parabolic barriers does not introduce additional param- 
eters with respect to the rectangular profile. 

The effect of an applied electric field £ is exemplified 
in Figure [TJ lower panel: the barriers are lowered and 
distorted, and eventually undergo a breakdown (regime 
which is out of the scope of the present paper). The 
energy profile changes to: U!(x) = Ui(x) — e£iX, where 
£i = Vi/di is the electric field within the i— th barrier, 
where the voltage drops by V%. 

The voltage induced change of shape determines signifi- 
cant changes in both thermal and tunnel currents. The 
thermal current density across the i— th barrier can be 
written as^ 




FIG. 1: Upper panel: sketch of the energy barrier profile 
used to describe the linear (small voltage) p c . The abscissa is 
the distance along the crystallographic c axis. Middle panel: 
modified energy profile chosen for the nonlinear transport. 
Lower panel: effect of an applied electric field (voltage) on 
the barriers. Note that the two barriers undergo different 
changes (uneven voltage drops) . 



energy, and is^ 
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where D is the diffusion coefficient, and n — const is 
the volume density of carriers outside the barrier. This 
expression can be solved analytically-^ to yield a closed 
form for J th ,i- 

To write down the expression of the tunnel current one 
needs the transparency of the i— th barrier. We in- 
corporate the effect of the electric field in analogy to the 
field effect. The final expression depends on the overall 
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1 otherwise. As- 



suming that the in-plane scattering affects the tunneling 
matrixes, and the transparency similarly, and that it holds 
in the nonlinear regime, the tunnel current is written as 
follows: 
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FIG. 2: Simulations of the dl/dV curves by means of the 
two-barrier model, using parameters typical for an overdoped 
BSCCO. Curves are calculated every 25 K for T < 300 K. 
With reference to Figure \T\ we used Ai = 2710 K, A2 = 
610 K. Panel (a): calculated curves arbitrarily shifted verti- 
cally. Panel (b): same curves, on absolute scales. Note the 
clear parabolic behavior and the emergence of the humps with 
lowering temperature. In real samples the crossing would be 
masked by the superconducting transition. 



where / is the Fermi function. 

To complete the model we need to stress that experi- 
ments are usually performed applying an external poten- 
tial V to the sample (typically a mesa crystal, composed 
of N cells). It is true that the total voltage distributes 
evenly across the N cells, so that it is customary to report 
the measurements as a function of the voltage per cell, 
v = V/N. However, within our model each cell is com- 
posed of two unequal barriers, so that v is the sum of the 
potential drops over the two barriers: v = V\ +V2, with 
V1/V2 variable with temperature or with v itself. The cal- 
culations must then be performed selfconsistently, with 
the requirement of the same total current through each 
barrier: J = J\ = J2, with Ji = Jth,i + Jtun,i (the cross 
section S is the same throughout the sample). 
Before presenting some of the results of the model, we 
recall that the model itself contains five adjustable pa- 
rameters: a, b, /?, Ai, A2, being p a b,n a measured quan- 
tity. None of the parameters are temperature dependent, 
while all can be doping dependent. Moreover, several 
constraints act to further reduce the freedom of choice 
of the parameters. First, a, b, (3 can be combined so that 



only two of them affect the shape of the calculated curves, 
the third acting as a mere scale factor. Second, the fits 
of p c at low voltage set a lower bound for Ai — A2. Once 
the set of parameters is fixed, one directly generates the 
full T-dependent family of curves dl/dV vs. V, and ob- 
viously p c (T). 

We would like to stress that the model sketched up to 
now is a pure "normal state" model: it takes into account 
thermal and tunnel transport across a series of two bar- 
riers between an otherwise conventional conductor, with 
the only specific introduction of the in-plane scattering 
effect on the transparency. In particular, this model does 
not introduce any specifically tailored DOS in the tunnel 
current, Eq. with implicit or explicit introduction 
of a gap opening at a certain temperature, and contains 
only the energies given by Ai, A2 and Ep . Moreover, 
the model does not concern the in-plane transport. In- 
stead, it takes p a \, from the experiments, assuming that 
it reflects the in-plane scattering (through e.g., the ratio 
of the parameters b/a). 

We now illustrate the predictions for the nonlinear dl/dV 
curves, presenting some simulations with sets of param- 
eters typical for overdoped samples and for underdoped 
samples. Following the common habit of presentation of 
the data, first we will present our simulations for dl/dV 
translating them arbitrarily on the vertical scale. By re- 
plotting the same curves on an absolute scale we will show 
that significant information is lost in the most common 
representation. 

In Figure we report the voltage dependence of the dif- 
ferential conductivity with parameters chosen to describe 
p c in an overdoped sample. In Figure^ we plot the data 
as it is customary, that is with some arbitrary vertical 
translation, while in Figure^) the same curves are replot- 
ted in an absolute scale. We first notice that the present 
model reproduces an otherwise unexplained feature of 
the experiments: the parabolic behavior of dl/dV in the 
normal state. In fact, this behavior is routinely observed 
in the experiments^^ but it is often ignored. Some- 
times the data are normalized to this "background",— 
attributing it to some instrumental feature or to self- 
heating (however, we note that the parabolic shape is 
observed also when self-heating is avoided by ns time- 
scale short-pulse measurements^). Clearly, such a nor- 
malization possibly removes some significant information 
from the data. 

A second striking result is the emergence of the well- 
known "humps" from the parabolic background. This is 
a very important result of our simulations: without any 
additional energy scale (such as a gap introduced in the 
DOS) we reproduce such humps only by virtue of the 
balance between the tunnel and thermal current across 
the barriers. Moreover, the simulations reproduce the 
nearly T-independence of the position of the hump, Vh, 
observed in many overdoped samples, and the increase of 
the hump height with decreasing temperature 

Changing parameters to those typical for underdoped 
samples several other features are revealed. In Figure 
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FIG. 3: As in Figure [2] for an underdoped sample. Curves 
are calculated every 25 K for T < 300 K. Barriers heights are 
Ai = 3270 K, A 2 = 1170 K. Panel (a): calculated curves arbi- 
trarily shifted vertically. Panel (b): same curves, on absolute 
scales. Note the crossing of all the curves at V x . 



[3] we report the simulations corresponding to underdop- 
ing. Again, several familiar experimental features emerge 
from the model. The humps still rise from the back- 
ground, with little temperature dependence of Vh, as in 
the previous simulations. Of great interest is the plot 
of the curves on absolute scale: we see that all the sim- 
ulated curves cross each other at a typical voltage V x . 
This feature is appreciable only in the absolute scale. It 
is noteworthy that this is exactly the experimental be- 
havior exhibited by underdoped mesas, as reported by 
many groups^iii 

Finally, note that in Fig.s 12131 the scales are compara- 
ble to those in typical experiments, thus suggesting that 
quantitative fitting is possible. 

The results of the simulations obviously call for precise 
measurements and direct comparison with the experi- 
ments, task accomplished in the following Sections. 



III. RESULTS AND DISCUSSION 

In this Section we compare our simulations with exper- 
imental data taken in small BSCCO mesas at Stockholm 
University^ Experimental details for data taken on sim- 
ilar mesas can be found in Reffl3l. 

One of the key points is the determination of the barrier 



heights and of the other parameters needed to reproduce 
both the linear (small voltage) p c (T) and the nonlinear 
differential conductance dl/dV curves. The parameters 
are determined as follows. Using the procedure illus- 
trated previously^ we fitted p c (T) in a large number 
of BSCCO single crystals at different doping S. The data 
were taken with an eight-terminal configuration^ or from 
the literature^ The resulting parameters exhibited very 
regular doping dependences, such as a linearly decreasing 
A2(S)J& By comparing the data taken in the mesa with 
the full body of the preexisting elaborations we identified 
the parameters to be used. The dl/dV curves were then 
calculated directly, without further adjustments. We re- 
mark again that we will compare the simulation and the 
data on an absolute scale, without any vertical transla- 
tions. 

On general grounds, we anticipate that the two-barrier 
model reproduces the linear and nonlinear behavior in 
the (supposedly) true normal state with great accuracy, 
while some quantitative discrepancy appear crossing the 
temperature T* where the pseudogap is supposed to be 
located. 

To investigate the pseudogap region we use the data 
taken on a severely underdoped sample with S <0.2135. 
Figure 0^ reports the comparison with the model. The 
calculations are performed using parameters appropriate 
to that doping. One can see that the shape of the calcu- 
lated curves is very close to the data, while the volt- 
age scale compared to the experimental data appears 
stretched. Interestingly, rescaling the voltage of the cal- 
culated curves by a factor 7 = 2 the agreement becomes 
impressive, as shown in Figure [4)3. We note that the 
qualitative features do not depend on the rescaling: the 
model reproduces the emergence and the evolution of the 
humps and, what it is particularly relevant, the crossing 
of the curves at a typical voltage. The voltage rescaling 
brings in quantitative agreement the model and the data. 
To appreciate the result we recall again the absolute ver- 
tical scale. 

The need for a voltage rescaling is a feature to be in- 
vestigated. The best approach is to compare the model 
to data in the slightly underdoped sample, where all the 
characteristic features are present: parabolic character- 
istics at high T, emergence of the humps, crossing of 
the curves in the "hump" region. Figure [5] reports the 
comparison of the model with the data in a sample with 
S =0.24. As it can be seen in Figured the parabolic be- 
havior is correctly reproduced, again without any volt- 
age rescaling. However, when the humps emerge and 
the crossing voltage appears, the model describes (nearly 
quantitatively) the data only after a voltage rescaling by 
a factor 7 « 1.5, similar to the factor needed in the 
severely underdoped sample. Noteworthy, the appear- 
ance of the scaling factor is connected to the apperance 
of the humps (and of the crossing voltage). We stress 
again that no modification of the DOS is included in the 
model: the evolution from parabola to humps appears 
only as a consequence of the balance between thermal 
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FIG. 4: Comparison of the model to the experimental data in 
an underdoped sample (S <0.2135). Panel (a): data and fits 
are reported on absolute scales. Barriers heights as in Figure 
131 Panel (b): the voltage scale for the fits has been divided 
by 7 =2. By the simple rescaling of the voltage scale the 
fits reproduce the data well. The qualitative features are all 
contained in the model. 



FIG. 5: Comparison of the model to the experimental data in 
a slightly underdoped sample (8 =0.24). Panel (a): data and 
fits are reported on absolute scales. Note that the parabola 
at high T is reproduced. The barriers heights are: Ai = 2880 
K, A2 = 780 K. Panel (b): the voltage scale for the fits has 
been divided by 7 —2. Rescaling of the voltage brings the fits 
in agreement with the data below the temperature where the 
humps appear, but the fit of the parabolic behavior is lost. 



and tunneling transport, and of the dynamically adjust- 
ing voltage drop between the two barriers. Moreover, the 
full dl I dV are fitted, in absolute voltage and conductance 
scale. This result is even more remarkable since only T- 
independent fitting parameters are involved. 

It can be concluded that several typical features of the 
data (in particular, the parabolic dl/dV and the emer- 
gence of the humps) are reproduced by the two-barrier 
calculations. When the humps appear, we find that the 
calculated dl/dVs require a scaling of the voltage by a 
factor 7 >1 to match quantitatively the data. By further 
lowering the temperature, the rescaled simulations track 
the position of the experimental features (humps, cross- 
ing), even if the experimental dl/dVs show a reduced 
dynamics with respect to the simulation (as expected in 
a pseudogap state). We believe that, while the humps 
are a feature inherent in the layered structure of BSCCO 
and, as such, captured by the simple two-barrier model, 
the additional physics that develops in the pseudogap 
state needs a specific incorporation into the model. The 
role of the voltage scaling factor is still unclear, but it 



is evidently connected to the electronic transformation 
at T*. To the best of our efforts, the rescaling cannot 
be achieved changing the other parameters in the tunnel 
and thermal current. We mention that, on purely nu- 
merical grounds, a rescaling of the electric charge e — > je 
produces the desired rescaling of the curves. However, 
before drawing strong conclusions on the charge of the 
carriers much more work would be needed. As an exam- 
ple, it should be considered that any "preformed pair" 
scenario would certainly affect the DOS and ultimately 
change the shape of the dl/dVs. We shall leave further 
work to future studies, but we emphasize again how a 
simple, purely "normal-state" model as the one here de- 
picted can reproduce many features of the experiments 
with surprising success. 
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IV. SUMMARY 

Summarizing, we have extended the two-barrier analy- 
sis of the c-axis transport for the normal state of BSCCO 
to the dl/dV curves, taken in mesas at different doping 
level. We have found that many of the peculiar features 
of such curves are reproduced by the two-barrier model, 
namely: (a) the parabolic dl/dV in the high-T region 
(above T*), (b) the emergence of the "humps" from this 
parabolic behavior and their nearly constant position (on 
the voltage scale) lowering the temperature, and (c) the 
crossing of the absolute dl/dV curves at a characteristic 
voltage V x . Our findings indicate that "conventional" 
transport across the energy barriers (tunneling and ther- 
mal activation) could be at the origin of most of the un- 
common features of the c-axis transport in BSCCO. In 
the normal state (above T*) simple tunneling and ther- 



mal activation reproduce the experimental dl/dV quan- 
titatively. Below T* quantitative discrepancies appear. 
A rescaling of the charge of the carriers and a reduction 
of the dynamic range bring the calculations in quantita- 
tive agreement with the experiments. The clarification of 
the normal state in terms of single electron tunneling and 
thermal activation could help to identify the possible sce- 
narios for the physics of the pseudogap state. The need 
for such modifications remains an open question that we 
intend to tackle in a future work. 
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